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Abstract 

^ . A simple method of the calculation of neutrino transition probabil- 

ities in vacuum in the general case of n massive neutrinos is presented. 
The method proposed fully utilizes the unitarity of the mixing matrix. 
Three-neutrino case for both neutrino mass hierarchies is considered 
in some details. Transitions in the case of the sterile neutrinos are 
^ also discussed. 

1 Introduction 

The observation of neutrino oscillations in atmospheric Super-Kamiokande 
P], solar SNO [2], reactor KamLAND [3] and other neutrino experiments jl] 
is one of the most important recent discovery in the particle physics. 

Neutrino oscillations are based on the assumption that states of flavor 
neutrinos i/ e ,i/^,u T and sterile neutrinos v Sl ,v S2 , ... are described by the co- 
herent superpositions of the states of neutrinos with definite masses (see ,for 
example, [5], ©, [7], [8], [9]) 



> = £^k>- (i) 
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Here n = 3 + n s (n s is the number of sterile neutrinos) , U is an unitary n x n 
mixing matrix, |z^) is left-handed state of neutrino with mass mj. The states 
\v a ) satisfy the condition 

(Va'Wa) = 8a>a- (2) 

If at t — the flavor neutrino v a is produced at the time t we have 

K)t = ^\vc/)Me~ iHot K), (3) 

where H is the free Hamiltonian. 
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Thus, the probability of the transition u a — > zv during time t is given by 
the expression 



iEitjT* 1 2 
ai I ' 



(4) 

It is obvious that the factor C/*j is the amplitude of the transition from 
initial flavor state into the state with definite mass the factor e~ lEit 
describes propagation in this state and the factor U a >i is the amplitude of the 
transition from the state with definite mass into the final state |zv). For the 
ultrarelativistic neutrino we have 
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E i=1 Jp' + Tr$~E + ^, (5) 

where E = p is the energy of neutrino at m, -> 0. 
The expression (JJJ) can be rewritten in the form 

P(v a -> vj) = | U a/i e- 2iA *U*J 2 . (6) 

i 

Here p is an arbitrary fixed index, 



A pi = (E l -E p )t^—^, Am 2 pi = m 2 - m 2 p . (7) 



Am^L 
AE 

We used the relation 

t = L, (8) 

where L is the source-detector distance^ 

From ([6]) it obvious that transitions between different neutrinos are pos- 
sible if two conditions are satisfied 

1. At least one neutrino mass-squared difference is different from zero. 

2. Neutrinos are mixed (U ^ 1). 



-"Xet us notice that the validity of the relation (JSJ) was confirmed in the high-accuracy 
recent OPERA measurement [10] 
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2 Standard expression for the transition prob- 
ability 

From ([6]) we obviously have 

P{v a -> v*) = \Ua'i\ 2 \U ai \ 2 + 2 ReJ^Ua'iKiU^U^e- 2 ^. (9) 

i i>k 

From the unitarity of the matrix U f/^ = <5 a / a ) we find 

X] \U aH \ 2 \U ai \ 2 + 2 Re^ U aH U* ai U*, k U ak = 5 a , a (10) 
j «>fc 

From and f llOp for the transition probability we obtain the following 
expression 

P[v a v*) = - 2 Re J] U ali U* ai U*, k U ak (l - e" 2iA -)- (11) 

«>fc 

Finally we obtain the following standard expression for the probability of 
neutrino transition in vacuum (see [8], [TT] ) 

P(z/ Q ->■ iv) = 5 Q / Q -4 ^ReUa'iU* iU^^Uak sin 2 A ki +2 } y ImUgfjU^U^^Ugk sin 2A fc j 

i>fc j>fc 

(12) 

In order to obtain z/ Q — > v a i transition probability we need to make in ffl2|) 
the change U ai —> U ai . Thus, we have 

P{v a ->■ *v) = 5 a / a -4 ^Re Ua'iU^U^^Uak sin 2 A ki -2 ^lmU a/i U* i U*, k U ak sin2A ki 

i>k i>k 

(13) 

It is obvious that in the case of a' = a the last terms of ( TT2"j) and ffT31) are 
equal to zero. We have 

P(y a ^v a ) = P{v a ^v a ). (14) 

This relation is a consequence of the CPT invariance. 

If CP invariance in the lepton sector holds, in this case U ai = ?7* f and 

P{y a v a t) = P{p a ~> v a >) a' 7^ a. (15) 
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From (I12p and f)13p for the CP asymmetry we have 

A C Ja = Pi"* -> Va>) ~ P("a *v) = 4 ^ A% a sin 2A fci . (16) 

i>k 

where 

4} a = Im U aH U* ai U* a , k U ak . (17) 
It obvious from (ITT)) that 

4J a = "^a- (18) 

In the case of the three- neutrino mixing the CP asymmetry is characterized 
by A 2 Ji a , A 3 J a and Af, a . We will show now that (see [12] ) 

In fact, taking into account the unitarity of the mixing matrix we have 

£4k = °- (20) 

From (TIHD and fl2QJ we find 

4 12 i 413 _n 421 ,423 _ n 431 , 4 32 _ n 

From these equations we easily find the relations ( Fl9|) . Thus, in the three- 
neutrino case the CP asymmetry has the form 

A c J a = 4 A 2 a ) a (sin 2A 12 + sin 2A 23 - sin 2A 13 ) (22) 

where 

A 13 = A 12 + A 23 . (23) 
For any a and b we have the relation 

■ 7 •/ i\ , . a + b . a . b . . 

sm fl + sino - sm(a + o) — 4 sin — - — sm — sm — . (24) 

Thus, the CP asymmetry in the three-neutrino case is given by the expression 
A c J a = 16 A 2 Q ) Q sin(A 12 + A 23 ) sin A 12 sin A 23 . (25) 
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3 Alternative way of the calculation of the 
transition probability 

We will present here a simple method of the calculation of the probability of 
the neutrino transition in vacuum. In expression for the vacuum transition 
probability, presented below, the unitarity of the mixing matrix will be fully 
utilized. In particular, the expression (125|) for the CP asymmetry will be 
obtained directly from the general expression for the transition probability 
without any additional calculations. 

Let us return back to the expression ([6]). Taking into account the unitarity 
of the mixing matrix we can rewrite this expression in the form 

P(y a _> Va ,) = \8 a , a + £ U^i (e~ 2iA - - 1) U*J 2 

i 

= \8 a , a -2i^2 U a >i U* ai e~ iA ^ sin A pi \ 2 (26) 

i 

It is obvious that the index i in (1261) runs over values % ^ p. 
From (1261) we have 

P(u a -)• v a ,) = 8 a > a - 4^ |i^| 2 sin 2 A pi 5 a > a + 4^ \ U a >i\ 2 \U a >i\ 2 sin 2 A pi 

i i 

+ 8 Re U aU U* m U* alk U ak e-^- A ^smA pi sinA pk (27) 

i>k 

Finally, we find the following general expression for v a — > v a i [p a — )■ v a i) 
transition probability 

P(u a -+v a >) = 5 a < a - 4^ \U ai \ 2 (5 a , a - |[/ Q /i| 2 )sin 2 A pi 

i 

+8 Re Ua'iU^U^^Uak cos(A pi - A pk ) sin A pi sin A pk 

i>k 

±8 Y Im U a'iU* ai U* a , k U ak sin(A pi - A pk ) sin A pi sin A pk (28) 

i>k 
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4 Three-neutrino oscillations 



4.1 General expressions for v\ — > vy {y\ — > Vy) transition 
probabilities 

In the case of the three-neutrino mixing there are two independent mass- 
squared differences. From analysis of neutrino oscillation data it follows that 
one mass-squared difference is much smaller than the other one. Correspond- 
ingly, two three-neutrino mass spectra are possible 

1. Normal hierarchy (NH) 

mi < ni2 < m,3, Am^ <C Am^. (29) 

2. Inverted hierarchy (IH)@ 



m 3 < mi < ni2, Am 2 2 <C |Am 2 3 |. (30) 

Let us denote two independent neutrino mass-squared differences Am| (so- 
lar) and Am 2 A (atmospheric). We have 

Am 2 12 = Am 2 s , m 2 23 = Am 2 A {NH) Am\ 2 = Am|, Am 2 13 = -Am 2 A (IH). 

(31) 

In the case of the NH it is natural to choose p = 2. From the general 
expression (I2"8"j) we have 

P NH (p\ V4) = 8 VI -AY,\Uii\ 2 {5 vl - \U V1 \ 2 ) sin 2 A 5 

i 

\Ui 3 \ 2 {Sm - \U l/3 \ 2 ) sin 2 A A - 8 Re U^U^U^Un cos(A A + A s ) sin A A sin A 5 

i 

T8 Im UnUyjyjUn sin(A A + A s ) sin A A sin A 5 (32) 

2 Notice that neutrino masses are labeled differently in the case of NH and IH. This 
allows to introduce the same mixing angles in both cases 
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In the case of the IH we choose p = 1. For the transition probability we 
obtain the following expression 

P IH $ V^) = 5 n - \Ui2\ 2 (S vl - \U, 2 \ 2 ) sin 2 A s 

i 

-4^|^ 3 | 2 (^ - |^ 3 | 2 )sin 2 A A -8 Re U V3 U^ 2 U l2 cos(A A + A s ) sin A A sin A s 

i 

±8 Im U V3 U; 3 U;, 2 U l2 sin(A A + A 5 ) sin A A sin A s (33) 

i>k 

The expressions ( 1321) and (|33|) differ by the change Uu — > U\ 2 and by the sign 
of the last term. Notice that for the CP asymmetry directly from (132]) and 
(|33|) we have 

A^f = -16 Im U^Uyj^Un sin(A A + A s ) sin A A sin A s (34) 
in the case of NH and 

= 16 Im U^U^U^Un sin(A A + A s ) sin A A sin A s (35) 



in the case of IH. In the standard parameterizations the 3x3 PMNS [13] .[14] 
mixing matrix U is characterized by three mixing angles and one CP phase 
and has the form 

C13C12 C13S12 si 3 e~ tS 

U = ( — C23S12 — s 23 c\ 2 Si 3 e lS C23C12 — s 23 S\ 2 Si 3 e lS C13S23 | . (36) 

■523 s 12 — C 23 Ci 2 Si 3 e lS — S23C12 — C 23 Si 2 Si2,e lS C13C23 

Here c\ 2 = cos 6*12, S12 = sin 6*12 etc. 
4.2 Leading approximation 

From analysis of the neutrino oscillation data it follows that two neutrino 
oscillation parameters are small: 

Am? 

— -4~3-l(T 2 , sin 2 9 13 ~ 2.4 • 1(T 2 (37) 
Am A 

In atmospheric region of the parameter ^ ( > 1) effects of neutrino 
oscillations are large. In the first, leading approximation we can neglect 
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small contributions of Am| and sin 2 # 13 into neutrino transition probabilities. 
From (I3"2"j) . (T3H]) and (1551) for the probability of (P M ) to survive (for both 
neutrino mass spectra) we obtain the following expression 



P NH & V^) ~ P^fe 5 ^ 1 - 4|f/ M3 | 2 (l - l^sl 2 ) sin 2 Ami 



4P 



1 -sin 2 2^3 sin 2 Ami — . ( 38 ) 

4P 



P(^ -> i/ T ) ~ 1 - P{p ll -> ^) ~ sin 2 2# 23 sin 2 Am^ — . (39) 



In the leading approximation we have Piy^ — >• z^e) — and 

L 

l 4£' 

Thus, in the atmospheric region predominantly two-neutrino z/,- oscil- 

lations take place. 

Let us consider now v e — > v e transition in the reactor Kamland region 
( ~ 1) • Neglecting contribution of sin 2 13 we have 

P NH {v e -> v e ) ~ P /H (z/ e ->■ P e ) ~ 1 - sin 2 2^2 sin 2 Am|-^. (40) 

4P 

For appearance probabilities we find 

P NH {v e -> P M ) ~ P /H (/> e ->• P M ) ~ sin 2 2# 12 cos 2 6 23 sin 2 Am|-^= (41) 

and 

P NH (D e -> P T ) ~ P m (^ e P T ) ~ sin 2 26> 12 sin 2 # 23 sin 2 Am|— (42) 

4P 

We have 

P{v e -»• P e ) = 1 - P(P e -> P M ) - P(P e -> P T ) (43) 

and 

P(p e ->■ i/ T ) 



tan 2 ^ ~1. (44) 



Thus, in the reactor Kamland region P e P M and v e v T oscillations take 
place. 

The expressions fl38l) and (1381 were used for analysis of the first Super- 
Kamiokande atmospheric data, K2K and MINOS accelerator data and data 
of the reactor KamLAND experiment. Now with improved accuracy of the 
neutrino oscillation experiments it is more common to perform more compli- 
cated three-neutrino analysis of the data. 
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4.2.1 v e — > v e survival probability 



From fl32|) and fl33|) we can easily obtain exact three-neutrino expressions 
for v e — > v e survival probabilities for both neutrino mass spectra. We have, 
corresp ondingly, 

P^(u e P e ) = 1 - 4 \U e3 \\l - \U e3 \ 2 ) sin 2 Am ^ L 



AE 



-4|t/ el | 2 (l-|^| 2 ) sin 2 ^^ 

4.C/ 



and 



|^3| 2 |^| 2 cos + Am|)L sin ^ sin (45) 

1 il 11 AE AE AE 



P IH (P e -+ p e ) = 1 - 4 |^ 3 | 2 (1 - l^esl 2 ) sin 2 A/// ' " 



4£ 



-4|£4 2 | 2 (l-|f/ e2 | 2 ) sin 2 ^ 



|9| , 9 (Ami + Ami) L AmlL AmlL , , 
-8 l^n^l 2 cos 1 ^ ^ sin sin (46) 

In the standard parameterization of the PMNS mixing matrix we have 
F Ns f - _^ D ) = 1 _ gin 2 w -2 Am^L 



— (cos 2 6*i3 sin 2 26> 12 + sin 2 26* 13 cos 6*i 2 ) sin 2 



2E 

AmlL 



2E 



9 „ 9 „ (Ara 2 j + Am 2 ,)I Am 2 A L AmlL , „. 

-2 sin 2 2#i 3 cos 2 6 12 cos - ^— ^— sin ^ sin — (47) 

4£< 4i? 4i? 



and 



pis- — ^ y\ = \ — sin 2 20n sin 2 



Am 2 L 

13""'" 



2£ 

Am 2 L 

13 Sill Zt7i2 "I" Sill Zt7i3 Sill t7 12 1 



— (cos 2 0i 3 sin 2 2019 + sin 2 20 13 sin 4 9 12 ) sin 2 



2£ 

9 „ 9 „ (Am 2 < + Am 2 ? )L AmlL AmlL , , 
-2 sin 2 2fl 13 sin 2 12 cos - ^ ^— sin ^- sin %-. 48 
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(-) (-) 

v n — >• v e vacuum 



4.2.2 — > v e (y^ — > is e ) transition probabilities 

From 022]) and we obtain the following expressions for 
transition probabilities: 

P NS K->%=4|[/ e3 | 2 |^ 3 | 2 sm 2 A A 
+4 |t/ e i| 2 |f/ Ml | 2 sin 2 A s 

-8 Re VaV^VlJJpx cos(A A + A s ) sin A A sin A s 

T8 Im UezU^U^U^ sin(A A + As) sin A A sin A s . (49) 

and 

P IS K = 4 \U e3 \ 2 \U^\ 2 sin 2 A A 

+4 |/7 e2 | 2 |f/ M2 | 2 sin 2 A s 

-8 Re U e3 U^U* 2 U^ cos(A A + A s ) sin A A sin A s 

±8 Im U e3 U; 3 U* e2 U^ sin(A A + As) sin A A sin A s . (50) 

Using the standard parameterization of the PMNS mixing matrix in the case 
of NH we have 

P NS K = sin 2 2^4, sm 2 

Ani 2 

+ (sin 2 29i2C 2 3 cl 3 + sin 2 2#i 3 Ci 2 s 23 + -f^c 2 2 coscT) sin 2 — — 

4i? 

9 x (Am\+AmV)L Am 2 A L Am%L 
+ (2 sin 2 2# 13 s 2 3 c 2 2 + K cos 5) cos ^ ^ ^— sin sin -^J- 

(Am 2 , + Aml)L Am\L AmlL , , 

=Fi^sin5 sin- — sin ^— sin . (51) 

AE 4E 4E v 7 

Here 

K = sin 2# 12 sin 2# 13 sin 2# 23 ci 3 . (52) 
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In the case of the inverted neutrino mass hiearchy we find the following 
expressions for the transition probabilities 

P IS K-^)= sin 2 28 13 s 2 23 sm 2 

4:11/ 

+ (sin 26 l i 2 c 1 3C 2 3 + sin 26 , i3S 12 s 2 3 — Ks 12 coso) sin — — — 

AE 

+ (2 sin 2 2^ 13 s 2 3 s 2 2 - AT cos 5) cos ^ ^ ^— sin sin 

(Am 2 A + Aml)L AmU AmU , , 

=FK sin 5 sin — sin — sin — . (53) 

AE AE AE y ' 

Formulas ( 1ST]) and f )53|) can be used for analysis of the data of T2K long 
baseline accelerator experiment in which matter effects are negligible. 



5 Transitions in the case of sterile neutrinos 

From existing data some indications in favor of neutrino oscillations driven 
by "large" (~ leV 2 ) neutrino mass-squared difference(s) were obtained. (see 
[15], [H]). These data (if correct) would proof existence of sterile neutrinos. 

From general formula (|28p we can easily obtain transition probabilities 
in such cases. In the framework of mixing of four massive neutrinos we will 
consider first the simplest 3+1 scheme in which the forth mass is separated 
from three close masses by a ~ leV gap. Let us choose p = 1. In the region 
of -| sensitive to large neutrino mass-squared difference ( A T E L ^ 1> Am 2 = 
Am 2 4 ) we have A^ ~ 0, % = 2, 3. From (j28p we find in this case 

P(u a Vzv) ~ 8 a , a - A\U a4 \ 2 {5 a , a - \U a , A \ 2 ) sin 2 A ^ L . (54) 
In the case of reactor antineutrinos we have 

P(u e ^u e ) = l-A \U ei \ 2 (l - \U ei \ 2 ) sin 2 (55) 



AE 



For v„ —7- v e transition (LSND) we find 



P(v» -> * e ) = 4 \U eA \ 2 \U^\ 2 ) sin 2 (56) 
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Let us consider now 3+2 scheme with 2 masses m 4 and m 5 separated from 
three close masses by a ~ leV gaps. Let us choose p = 1. In the region of -| 
sensitive to Am 2 4 and Am 2 5 we have Aij ~ 0, i — 2, 3. From ( 128]) we find 
in this case 

PR = 4|[/ e4 | 2 |^ 4 | 2 ) sin 2 + 4|^ e5 | 2 |[/ M5 | 2 ) sin 2 

+8 Re U^IT^ cos(^| ^) sm -JL_ sm 

TT „ .Am 2 t:L Am?,L AmU Am?,! , , 
±8 Im U^U*^ sm(^| ^) sm sm (57) 

FotVq — >Vq survival probability we find the following expression 

pRVzt) = l-^^ra-i^nsin^^-^^m-i^i^sm^^ 

|9| .„ ,Ara?,L Am?,L, Am?,L Am?.L , , 
+ 8^ 2 |£U 2 cos(^f i^) sin ^ sin ^ ^ 

I am thankful to A. Olshevskiy and C. Giunti for useful discussion. 
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